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Abstract 

Burgers vortices are stationary solutions of the three-dimensional Navier-Stokes 
equations in the presence of a background straining flow. These solutions are given 
by explicit formulas only when the strain is axisymmetric. In this paper we consider 
a weakly asymmetric strain and prove in that case that non-axisynnnetric vortices 
exist for all values of the Reynolds number. In the limit of large Reynolds numbers, 
we recover the asymptotic results of Moffatt, Kida & Ohkitani m • We also show 
that the asymmetric vortices are stable with respect to localized two-dimensional 
perturbations. 


1 Introduction 

Localized structures such as vortex sheets or tubes play a prominent role in the dissipation 
of energy in three-dimensional turbulent flows. It is believed that these dissipative struc¬ 
tures take place due to the interplay of two basic mechanisms: amplification of vorticity 
due to stretching, and diffusion through the action of viscosity EQ- A typical example 
that exhibits both features is the familiar Burgers vortex jjlj, an explicit solution of the 
three-dimensional Navier-Stokes equations in the presence of an axisymmetric background 
straining flow. In real flows, however, the local strain has no reason of being axisymmet¬ 
ric, and as a matter of fact the vortex tubes observed in numerical simulations usually 
exhibit a truly elliptical core region. It is therefore important to study the analogue of 
the Burgers vortex when the straining flow is asymmetric, although no explicit expression 
is available in that case. 

Using a double series expansion, Robinson and Saffman H3 formally established the 
existence of an asymmetric vortex for small values of the Reynolds number R and of the 
asymmetry parameter A. This solution was also studied numerically for larger A (up to 
3/4) and R (up to 100). On the other hand, an asymptotic expansion for large Reynolds 
numbers was performed by Moffatt, Kida and Ohkitani mi, see also 0. Their results 
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indicate that an equilibrium stretched vortex should exist for all values of A G (0,1) and 
R > 0 such that \/R -C 1. Interesting features of these solutions, such as the shape of 
isovorticity contours and the spatial distribution of energy dissipation, were also studied 
in detail. Finally, the stability of symmetric or non-symmetric vortices is an important 
issue which has attracted a lot of attention in recent years. Roughly speaking, the stability 
with respect to two-dimensional perturbations (i.e., perturbations which are independent 
of the axial coordinate) is well understood [T7] [T3, EJ [7j, but only partial results have 
been obtained in the general case where arbitrary three-dimensional perturbations are 
allowed nail El IH El 

In this paper, we prove (rigorously) that non-axisymmetric Burgers vortices exist 
for all values of the Reynolds number, provided the asymmetry parameter is sufficiently 
small. In particular, taking the limit R —> oo, we recover exactly the asymptotic results of 
Moffatt, Kida and Ohkitani mi- We also show that these vortices are stable with respect 
to spatially localized two-dimensional perturbations. Existence for larger values of the 
asymmetry parameter and stability with respect to three-dimensional perturbations are 
difficult questions, which have been solved so far for small Reynolds numbers only [ 8 J. 

We now describe our results in more detail. We consider an incompressible viscous 
fluid Elling the whole space R 3 , and we suppose that the velocity field is a two-dimensional 
perturbation of a linear straining flow, namely 

/ 71^1 \ /u 1 (xi,x 2 ,t) \ 

U(xi,x 2 ,x 3 ,t) = 72 X 2 + u 2 (x 1 ,x 2 ,t) , 

\ 73^3 / \ 0 / 

where 7 i, 72,73 are reals constants satisfying 71 + 72 + 73 = 0. Throughout this paper we 
assume that 

7i = -|(1 + A), 72 = -|(l-A), 73 = 7, (1) 

for some 7 > 0 and some A G [0,1). Thus 73 is the only positive principal rate of strain, 
and the straining flow is axisymmetric if and only if A = 0. The case of a biaxial strain 
(A > 1), which is also important for applications in turbulence, will not be considered 
here. 

The vorticity fl = V x U is aligned with 
horizontal variable, namely 

( 0 \ 

n(x 1 ,x 2 ,x 3 ,t) = 0 

\u(xi,x 2: t) J 

The evolution equation for oj reads 

d t oj + (ui + 7 iaq)di oj + (u 2 + r j 2 x 2 )d 2 oj = uAuj + 73 a; , ( 2 ) 

where v > 0 is the kinematic viscosity of the fluid. Since d\U\+d 2 u 2 = 0 and d\u 2 — d 2 u\ = 
cu, the rotational velocity u = {u\,u 2 ) can be recovered from u via the two-dimensional 
Biot-Savart law 

U 0M) = 7 T [ 7 -— ^wu(y,t)dy, x = ( 27 , x 2 ) G R 2 , (3) 

Z 7 T ,/ R 2 \x — y | 


the vertical axis and depends only on the 


where oj = d\U 2 — d 2 U\ . 
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where x L = (—x 2 ,x 1 ) and \x\ 2 = x\ + x\. 

In the axisymmetric case A = 0, equation © has a family of explicit time-independent 
solutions: 


a/ = 


<5 2 


where T e R, 5 = (V/ 7 ) 1 / 2 , and 


G(X) = 


~| X | 2 / 4 ) = 


u r = 


1 x 1 
27r \x\ 


r 

~s 


r G 


(X\ 

KSJ ’ 


(4) 

-W 2 /4A ) 

x e R 2 . 

(5) 


These are the well-known (axisymmetric) Burgers vortices. The family is indexed by the 
parameter T = f R2 u; r da:, which represents the circulation of u r at infinity. The Reynolds 
number associated to the Burgers vortex with circulation T can be defined m as 


R = — . 


The aim of this paper is to study the analogue of the Burgers vortices when the 
straining flow is not axisymmetric. The expressions of these asymmetric vortices will be 
greatly simplified if we use the natural lengthscale 5 = (v/ 7) 1 ^ 2 and timescale r = I /7 
defined by the viscosity and the strain. We thus replace the variables x, t and the functions 
u, u with the dimensionless quantities 

x ~ t ru 

X — — , t = - , U = —— , U = TLO . 

or 0 

Dropping the tildes for convenience, we see that the new functions u>, u satisfy O with 
7 = v — 1 , namely 

dtUJ T u • Vo; — Huj T , (6) 

where 

£ = A + |x-V + l, M = \{xidi - x 2 d 2 ) . (7) 

It is easily verified that LG = 0 and v G ■ VG = 0, hence in the symmetric case A = 0 the 

Burgers vortex u a = aG, u“ = av c; is indeed a stationary solution of © for any a G R. 

The Reynolds number associated to this flow is simply R = |a|. 

To formulate our results, we introduce appropriate function spaces. Let X be the 
(real) Hilbert space 


X = |u> € L 2 (R 2 ) 
equipped with the scalar product 

(wi,w 2 ) x ■ 

We also define the subspace Y = 
natural scalar product 


G~ 1/2 w e L 2 ( R 2 ), 




( 8 ) 


in : 2 


G(x) 1 w 1 (x)w 2 (x) dx 


{w E X | diW G X for i = 


1 , 2 } equipped with the 


(wi,w 2 )y = J G(x) 1 ywi(x)w 2 (x) + Xw\(x) ■ Xw 2 (x)J dx . 

Except for the zero mean condition, the space A" is just a weighted L 2 space (with Gaussian 
weight) and Y is the corresponding Sobolev space. We can now state our Erst result: 
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Theorem 1.1 There exist A 0 > 0 and K 0 > 0 such that, for all A G [0, A 0 ] and all a G R, 
equation (GJ) has a unique stationary solution u; Q,A ,u“ ,A such that ||u;“’ A — ccG||y < K 0 . 
Moreover, ||u; Q ’ A — aG||y < A' 0 (A/Ao)|a|/(l+|a|). 

This theorem shows that, if the asymmetry parameter A is sufficiently small, equa¬ 
tion © has a family of equilibria a; a,A indexed by the circulation number a = f r2 u) a,x dx. 
The solution u;" ,A is locally unique, has a Gaussian decay at infinity, and converges to 
the Burgers vortex aG as A —> 0. Further properties of these asymmetric vortices will be 
established in Section El For instance, uj a,x (x) is a smooth function of x G R 2 , a G R, 
and A G [0, A 0 ]. Moreover, LU a,x > 0 if a > 0, u> a,x < 0 if a < 0, and LU a,x = 0 if a — 0. 

Theorem II. II will be proved by a classical perturbation argument. The only remarkable 
point is that this argument can be applied uniformly for all a G R. In particular, for 
fixed A, we can investigate the limit of large Reynolds numbers R = |a| —> oo. The 
last inequality in the theorem asserts that a~ x uj a,x = G + 0(X/R) as R —> oo, hence we 
recover the observation by Moffatt, Kida and Ohkitani nn that the asymptotic profile is 
always the Gaussian G, even in the asymmetric case A > 0. (Note that, unlike in |11 j . we 
do not need to assume that the asymptotic profile is radially symmetric.) Moreover, the 
deviation from the limiting profile is proportional to X/R at leading order, as established 
in HD- A rigorous expansion up to second order in A and R 1 will be performed in 
Section 0 see Ea. flHHl) below. 

Our second result shows that the asymmetric Burgers vortex is asymptotically stable 
stable with respect to perturbations in A". As in Theorem 11.11 this property holds uni¬ 
formly for all a G R and A G [0, Ai], for some Ai > 0 (possibly smaller than Ao). Remark 
that there is no loss of generality in assuming that the perturbations have zero mean, 
because if J Cj dx = /3 ^ 0 then uj a,x + u) is a zero mean perturbation of the (modified) 
vortex cu“ +/3,A . 

Theorem 1.2 Given any /j G (0,1/2), there exist X± > 0 and e > 0 such that, for all 
A G [0, Ai] and all a G R, the following holds. For all initial data uj 0 with \\u> 0 —uj a ’ x \\x < £, 
equation WH has a unique global solution u(x,t ) such that u — tv a,x G G°([0,+oo),X). 
Moreover, 

||u;(-, t) — u;"’ a ||a' < !|a>o — u a,x \\x , for all t > 0. (9) 

The proof of Theorem 1 1.1 1 is based on the following ideas. Given A G (0,1) and a G R, 
we look for stationary solutions of (ED of the form u = aG + w, u = «v G + v, where 
f R2 w dx = 0 and v is the velocity field obtained from w via the Biot-Savart law (J2D- The 
equation for w reads: 

a(v G ■ Vw + v • VG) + v ■ Viv = Cw + XM.(aG + w) . (10) 

Let A be the integro-differential operator defined by 

Aw = v G -Vw + v-VG . (11) 

It is shown in j7] that, for any a G R, the spectrum of C — aA acting on X is contained 
in the half-plane {z G C | Re (-) < -5>- In particular C — aA is invertible, and m can 
be rewritten as 

w = (£ — a:A) _1 ■ Vw — XA4(aG + w) j . (12) 
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In Section [21 we show that (£ — aA)^ 1 ^ is a bounded operator in Y whose norm is 
uniformly bounded for all a G R. In Section |U we prove that the function w a E Y 
defined by 

w a = —a(£ — aA )~ l XiG (13) 

is uniformly bounded in Y for all a E R, and converges to some limit as |a| —> cx). 
After these preliminaries, a standard contraction argument allows to prove that m has 
a unique solution w a,x (in an appropriate ball in Y) if A is sufficiently small, and that 
w a,x = Xw a + 0( A 2 ). This is done in Section 4, where additional properties of the 
asymmetric Burgers vortex w a,x = aG + w a,x are also established. Theorem 11.21 is proved 
in Section [21 by an energy estimate, using the observation that the linearization of © at 
the vortex a;“ ,A is a small perturbation of the linear equation d t uj — (£ — aA)u> if A is 
small. 


2 Linearization in the symmetric case 

If £, A are the linear operators defined in ©, (HU), we know from [7] that £ — aA is 

invertible in A" (with bounded inverse) for all a G R. In this section we use the methods 

of [Zj to establish the following result: 

Proposition 2.1 There exist positive constants K \, K 2 such that, for all a G R, the 
following inequalities hold: 

||(£ — ccA) _ 1 w||y < AAlMIx , for all w G X , (14) 

|| (£ — aA)~ l M.w\\Y < K 2 \\w\\ y , for allw EY . (15) 

Proof: The properties of the linear operator £ acting on X are easy to establish, because 
this operator is conjugated to the Hamiltonian of the harmonic oscillator in R 2 . Let H be 
the closed subspace of L 2 (R 2 ) defined by H = {/ G L 2 (R 2 ) | f R2 G 1//2 fdx = 0}. Consider 
the linear operator L : D(L) —> H defined by 

D{L) = {f E H \ Af G L 2 (R 2 ), \x\ 2 f G L 2 (R 2 )} , 

L = G~y\-C)G^ = -A + LI - 1 

lb 2 

As is well-known, L is self-adjoint in H with spectrum u(L) = {n/2 \ n = 1, 2, 3,...}. In 
particular, L >1/2. It follows that there exists C\ > 0 such that, for all / G H, 

\\L~ 1/2 f\\ L 2 + |||x|L ~ 1/2 /|| L 2 + ||VL ' 1/2 /|| L 2 < C ' 1 ||/|| L 2 . (16) 

Indeed, setting g = L -1 / 2 /, we have 

11 / 111 = = \\L 1,2 g\\h = (9-Lg) L , 

= \S^ v ^ 2+ ^ g2 -\ g2 ) ix+ \ (g ' Lg)L ' (i7) 
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where we have used the fact that (g,Lg ) L 2 > l||g||| 2 . The bounds (Hfll) also show that 
L -1 / 2 is a compact operator in H. 

Another useful estimate can be obtained from (ED by a duality argument. For j — 1 or 
2, let Hj — {/ E H f R2 XjG 1 / 2 fdx = 0}. Then L~ x l 2 Xj extends to a bounded operator 
from Hj into H, and for all / G Hj we have 

\\L- l ' 2 Xjf\\ L 2 < c x ||/|| l2 . (18) 

Indeed, by density, it is sufficient to prove ED for f £ Hj fl 5(R 2 ), where S denotes the 
Schwartz space of test functions. In that case Xjf E H and for all p E H we have 

\(ip,L- 1 / 2 Xjf) L 2\ = \(xjL- 1 / 2 p,f) L 2\ < \\xjL~ 1 / 2 ip\\ L 2\\f\\ L 2 < Ci|MMI/IIl 2 , 

which proves ED- 

We now return to the operator C : D(C) —> X defined by D(C) = G 1//2 .D(L), X = 
G l / 2 H, and — C = G 1 / 2 LG ^ 1 / 2 . By construction, C is selfadjoint in A", — C > 1/2, and 
(—£) -1 / 2 is a compact operator in X. Using (fTcTl) . (fTRli we easily obtain the following 
additional properties: 

Lemma 2.2 

i) (—£) -1 / 2 is a bounded operator from X into Y; 

ii) (—£) -1 / 2 A4 extends to a bounded operator from. Y into X. 

Proof of Lemma HOI 

i) Let w E X. Since (-£)“ 1/2 = G l ' 2 L~ l / 2 G ~ 1 / 2 , we have 

U-C)~ 1/2 w\\ x = \\G~ 1 / 2 (-jC)- 1 / 2 w\\ L 2 = \\L- 1 / 2 G- 1 / 2 w\\ L 2 

< Ci\\G~ l ^ 2 w \\ L 2 = CMx . 

Similarly, since VG 1//2 = — |G 1 / 2 , we fold 

\\X(~C)- 1 / 2 w\\ x = \\G~ 1 / 2 XG 1 / 2 L- 1 / 2 G~ 1 / 2 w\\ L 2 

< \\XL- 1 / 2 G~ 1 / 2 w\\ L 2 + i|| \x\L- l / 2 G~ l/ 2 w\\ L 2 

< Ci\\G~ 1 ^ 2 w\\l 2 = C 1 \\w\\ x . 

ii) For any w E Y, we have ||(— C)~ l ^ 2 Mw\\x = \\\L~ l l 2 G~ l l 2 {;xidi— X 2 d 2 )w\\ L 2 . Now 
G~^ 2 djW E Hj for j = 1,2, hence by (fTHl) 

||L _1 / 2 XjG _1 ^ 2 (? J w||i,2 < CiWG^^djivWtf < UiHwllr . 

We conclude that ||(—£) _1 / 2 A1w||.y < Ci||iu||y. □ 

Finally, we consider the operator £ = (—jC) -1 / 2 A(—£) _1 / 2 , where A is defined by ED- 
The following properties of £ will be useful: 

Lemma 2.3 

i) The operator £ is compact in X; 

ii) The operator £ is skew-symmetric in X; 

Hi) For any a E R. the operator 1 + a£ is invertible in X and 

||(1 + a£) _1 in||jY < 11in11 „y , for all w E X . (19) 
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Proof of Lemma 12.31 

i) Since (—£) _1 / 2 is compact in X and bounded from X into Y, it suffices to show that 
A : Y —> X is bounded. If w G Y, then 

||v G • Vw||x < j|v G |UMMHU < C||w||y . 


On the other hand, by Holder’s inequality, we have for all p G [1,2]: 

MU’ = < IIG ,/2 ll i * = dMI* . ( 20 ) 

If v denotes the velocity field obtained from w via the Biot-Savart law ©, the Hardy- 
Littlewood-Sobolev inequality H3 implies that v G L 9 (R 2 ) for all q G (2, 00 ), and the 
following bounds hold: 

MU'? < Cjw\\ LP , where 1 <p < 2 and - =-. ( 21 ) 

q p 2 


Thus 


||v VG||x = ||G“ 1 / 2 v- VG || L 2 < IMMM-^VGIM < CIMIU 4/3 < C\\w\\x . 

We conclude that ||Aw;||x < C'||w|M 

ii) It is shown in ([ 7 ], Lemma 4.8) that (wi, Aw 2 )x + (Awi,w 2 )x — 0 for all w±,w 2 G Y. 
Since (— C)^ 1 ^ 2 is symmetric in X, it follows that (irq, Hw 2 )x + (Swi,w 2 )x = 0 for all 
w u w 2 G X. 

iii) By the analytic Fredholm theorem we know that 1 + aE is invertible in X for 
all a G R except perhaps on a discrete set (with no limit point) where the meromorphic 
map ot 1 —> (1 T ctE ) -1 has poles. But whenever 1 + aX is invertible we have by ii): 

|Hlx = (w,(l+a£)w)x < |M||.y|| (l+aE)ic||x for all w G X , 

hence ||(1 + aE) _ 1 i(;||jY < ||w||x- This implies that 1 + aE is invertible for all a G R and 
that m holds. □ 

Equipped with these lemmas, it is now straightforward to conclude the proof of Propo¬ 
sition Io For any a G R, the formula 

(C-aA)- 1 = —(—£) -1/2 (l + aE ) _1 (—£)~ 1/2 

shows that C — aA is invertible in X for all a G R. The bounds m and m are then 
direct consequences of this identity and Lemmas [T21CZ3 □ 


3 Large Reynolds number asymptotics 

The main goal of this section is to prove that the function w a defined by m is uniformly 
bounded in the space Y for all a G R. From (THU) we expect that w a —> A^AiG as 
| a | —> 00 , but it is not clear a priori that this limit makes sense because A is not an 
invertible operator. Our first result shows that AiG is indeed in the range of A (we recall 
that 5(R 2 ) is the Schwartz space of test functions): 
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Proposition 3.1 There exists G 5(R 2 ) such that Aw^ = MG. 

Proof: By (0), © we have MG = |(xi(9 1 — x 2 d 2 )G = — \{x\ — x\)G. Using polar 
coordinates in R 2 , we thus find 

(MG)(r cosh, r sin 6 *) = — ^r 2 g(r) cos(20) , ( 22 ) 

where g(r) = G{r cos 6,r sin 9) = ( 47 r) _ 1 e _r2 / 4 . As was observed in jTTJ El Ij, the 
operator A is invariant under rotations in the plane, and is therefore block-diagonal in 
the Fourier basis { e in0 } n£ z■ We make the following Ansatz: 

Woo(r cos9,r sin 6 *) = o;(r) sin(26 ) ) , (23) 

where u : R + —> R has to be determined. The velocity field associated to reads 

2 

Voo = -Q(r) cos(26 ) )e r — 0 (r) sin(26 ) )e e , 
r 

where e r is the unit vector in the radial direction, e# = e^r, and where : R + —> R is the 
solution of the ordinary differential equation 

— -(rilTr))' + 4rfl( r ) = UJ ( r ) > r > 0 , (24) 

which satisfies the boundary conditions 0(0) = f2(+oo) = 0. Using these expressions, we 
find 

(Aw^)^ cos 9, r sin 9) = (v G ■ Vwoo + • VG)(r cos #, r sin 0) 

= cos(26*)(2^(r)a;(r) — g(r)Q(r)) , (25) 

where ip{r ) = (27rr 2 )^ 1 (l — e~ r If we compare (E2D and (ED we obtain the solution 

v(r) = h(r) (fi(r) - , (26) 

where h = g/(2cp), i.e. h(r) = (r 2 /4)(e r2//4 — 1) _1 . Inserting (©Hi into (1231) . we see that U 
should satisfy the ordinary differential equation 

-i(rfi'(r))'+(^-M0)!2M = -Afh, r>0 . (27) 

together with the boundary conditions 0 ( 0 ) = Q(+oo) = 0 . 

Remark 3.2 Equation was derived and studied numerically by Moffatt, Kida and 
Ohkitani, see Eq.(2.25) in J77J/. The notation used in m is f = VL and g = —h. 

To solve ED, we first consider the associated homogeneous equation 

--(rO'(r))'+ (l^( r ))n(r) = 

8 


0 , r > 0 . 


( 28 ) 



Setting r = and Vt(r) = F(±log(r)), this equation is transformed into 


- F"(t) + (4-H±(t))F{t) = 0 , te R, (29) 


where H±(t ) = e ± 2 *h(e ±i ). In particular, H±(t) decays rapidly to zero as t —> +cx). 
Applying Theorem 3.8.1 in ; 3], we deduce that © has a unique solution F±(t) such that 


lim e 2t 

£—»+oo 





We now define tp+pr) — F + (log(r)) and "0-( r ) — F_(— log(r)). By construction, 
ip + ,'ip- are the only solutions of (I2%1) such that 

ip+(r) ~ — as r —> +cx) , -0_(r) ~ r 2 as r —■> 0 . 

We observe that the “potential” 4/r 2 — h(r) in is strictly positive, because 

inf(— h(r)) = inff---—) > 0 . 

r>o\r 2 J zxAz e z — 1 / 

By the Maximum Principle ra. it follows that ip'_(r) > 0 and ip' + {r) < 0 for all r > 0. In 
particular, ip + and ip- arc linearly independent, hence there exists Wq > 0 such that 

W(r) = ip + (r)ip'_(r) — ip' + {r)^)-{r) = — , r > 0 . 

Moreover, 

, / x n , , v W 2 

ip+(r) ~ —- as r —> (J , W-\ r ) ~ — i — as r —> +oo . 


Using these notations and the “variation of constants” formula, we obtain the following 
expression for the solution of P71>: 


O(r) = —ip + (r) 


,r ,3 


4 w ( 


r i/)-(z)h(z) d z — ip-{r) 


z 3 


4 W( 


7p + (z)h(z) dz , r > 0 . (30) 


It is clear that : R + —» R is a smooth function satisfying 


U(r) 


Q + r 2 

Q_r 


-2 


as r —> 0 , 
as r —> +cx) , 


where = 




4ta f 


ip±(z)h(z) dz . 


(The values ff + ~ —0.38 and ~ —17.5 were found numerically in jTTj.) Similar 
estimates hold for all derivatives. Going back to (I2H1) , we see that u : R+ —> R is smooth 
and rapidly decreasing at infinity. Moreover u;(0) = 0, t o(r) < 0 for all r > 0, and it 
is easy to verify that the Taylor expansion of uj{r) at r = 0 contains even powers of r 
only. Thus the function the function : R 2 —> R defined by (1231) is smooth, rapidly 
decreasing at infinity, and satisfies A = M.G by construction. □ 

As a consequence we can prove that the function w a defined in m is uniformly 
bounded in Y: 
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Corollary 3.3 There exists K 3 > 0 such that ||w a ||y < A' 3 |a|/(l+|a|) for all a G R. 

Proof: From (THfl) and (jT51) we know that ||wa,||y < K 2 \a\ ||G||y. On the other hand, in 
view of Proposition 13. 11 we have (£ — aA)w a = —aM.G = (£ — aA)w oc — £w 00 , hence 

w a = Woo - (C - aA)~ l Cwoo • (31) 

Using (HD) , we infer that ||w a ||y < ||woo||r + i£i||£woo||x- Combining both estimates we 
obtain the desired result. □ 

A more detailed analysis reveals that Woo is indeed the limit of w a as |a| —> oo: 

Proposition 3.4 There exists K A > 0 such that \\w a — Woolly < AT 4 /( l+|a|) for all 
ol G R. 

Proof: By (I3TT) and (1141) . we have \\w a — Woo||y < Ailj/iwoollx- On the other hand, 
proceeding exactly as in the proof of Proposition 13.11 it is straightforward to show that 
there exists z oo G 5(R 2 ) such that A^oo = Cwoo (the details are left to the reader). Then, 
using m, we find for all a^O 

W Q = Woo + ~Zoc ~ -(£ - aA) _1 £^oo , 
a a 

hence ||w Q — Woo||y < |o| _ 1 (ll z oo||y + A'xH^oollx)- Combining both estimates we obtain 
the desired result. □ 


4 Existence of asymmetric vortices 

As is explained in the introduction, we shall prove the existence of a stationary solution 
of (P) by solving Eq. (HD, namely 

w = A w a + (£ — aA ) -1 ■ Vw — \M.w ^ , (32) 

where w a is defined in m and (as usual) v denotes the velocity field obtained from w 
by the Biot-Savart law ©• To bound the nonlinear term in (ED, we use the following 
bilinear estimate: 

Lemma 4.1 There exists K 5 > 0 such that, if w,w G Y and if v is the velocity field 
obtained from w by the Biot-Savart law, then v • Vw G X and 

||v ■ Vw|| X < AA||w||y||w||y . 

Proof: We first observe that w G L P (R 2 ) for all p G [1, +oo), and that ||w||lp < Cp||w||y. 
Indeed this is proved in ED if p G [1,2], and if p > 2 this follows from the embeddings 
Y FA(R 2 ) > L p { R 2 ). Next, by the Calderon-Zygmund inequality the velocity 
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field v associated to w satisfies || Vv|| iP < Cp||w||iP for 1 < p < oo. Using the Gagliardo- 
Nirenberg inequality H2, we deduce that v G L°°(R 2 ) and that 

IMU« < C-llVvIl^HvHi? < C|H|^ 2 |H |^ 3 2 /2 < C\\w\\ Y , 

where we also used the Hardy-Littlewood-Sobolev inequality (THT) with q — 6 , p — 3/2. 
We conclude that v ■ Xw G X and 

||v • Vw||x < ||v|| L oc||Vw||x < G||u;||y||w||y , 

which is the desired result. □ 

We can now prove the main result of this section: 


Proposition 4.2 Choose Ao G (0,1) such that 

2/l 2 Ao < 1 , and I 6 A/A 3 A 5 A 0 < 1 . (33) 

Then, for all A G [0, A 0 ] and all a G R, equation FM) has a unique solution w a ’ x G Y 
such that ||tc a,A ||y < (4/1/A/) W This solution depends smoothly on the parameters 
(cc, A) G R x [0, Aq] and satisfies ||tc“ ,A ||y < 4A/A|cx|/(l+|cx|). 


Proof: Fix A G [0, Aq], cc G R, and choose r > 0 such that 


4A ? A * ^ 


where Kf = 


KM 

1 I Oi I 


(34) 


Let F a,x : Y — > Y be the quadratic map defined by 

F a,x (w ) = A w a + (£ — cxA)' 1 • Xw — AMw'j , w G Y . 
Using Proposition 12.11 Corollary 13.31 and Lemma 14. 11 we find for all w G Y: 

||A“’V)||y < Kf\ + KJiMWy + K2X\\w\\ Y ■ 

Similarly, for all w , w G Y: 

||F“' a (to) - F“' J (ro)|| r < (k^sIWwWy + ||ro||y) + K 2 X\ ||to - ro||y . 


(35) 


(36) 


Let B r = {w G Y | ||u>||y < r}. It follows from El, El, and El that F q,a maps B r 
into itself, because A/* A + A/A/r 2 + A/Ar < r/4 + r/4 + r/2 = r. Similarly (13611 implies 
for all w, w G B r : 

||A"’ A (tc) — A a ’ A (-u))||y < k\\w — w\\y , where k — \ + K 2 X 0 < 1 . 

By the contraction mapping theorem, F a,x has thus a unique fixed point in B r , which we 
denote by w a,x . It remains to show that w a,x is a smooth function of (ex, A) G R x [0, Ao]. 
But this is a direct consequence of the implicit function theorem, because the map F a,x 
depends smoothly on (cx, A) and the differential 


D w F a,x (w ) = w i—► (£ — cxA) 1 • Xw + v ■ Xw — AM.w'j 
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satisfies \\D w F a,x (w)\\ < k for all w G B r and all (ck, A) G Rx [0, A 0 ]. Thus 1 — D w F a,x (w) 
is invertible at w — w a,x , and the desired conclusion follows from the implicit function 
theorem. □ 

Theorem o is an immediate consequence of Proposition 14.21 we just set ui a,x = 
aG + w a ’ x , u" ,A = av G + v" ,A (where v“ ,A is the velocity field obtained from w a,x by the 
Biot-Savart law), and K 0 = (AKiK^) 1 . By construction, u a,x is a stationary solution of 
© satisfying the conclusions of the theorem. □ 

In the rest of this section, we establish a few additional properties of the asymmetric 
vortex uj a,x \ 


1) Expansion in A. There exists K 6 > 0 such that, for all (a , A) G R x [0, Ao]: 

||w“’ A - aG - Xw a 11 y < As ^ A , 2 . (37) 

1 + |a| 

Indeed, using the notations of Proposition 14.21 we have w a,x G B r with r = 4A'^A. As 
w a,x is a solution of 432), we obtain 

|| w a ' x - Xw a \\ Y < KJQr 2 + K 2 Xr < K ^ a j X . 

1 + a 


2) Large Reynolds number asymptotics. Combining (TT7T) and Proposition 13.41 we find for 

|a| > 1: 


—u> a,x -G- - Woo 
a a 



k 6 a 2 k a x 

\a\ \a \ 2 


(38) 


In agreement with mi, we see that the leading correction to the Gaussian profile G is 
(A/a)woo, and that the higher order corrections are proportional to A 2 /R and A /R 2 , where 
R — |o;| is the Reynolds number. 

3) Small Reynolds number asymptotics. Since CAA.G = C~ l M.G = —A4G, it follows from 
m w a = aM.G+a 2 (C—aA)~ 1 AJ\4G. Replacing into (1371) and using Proposition ^. 11 
we obtain for |a| < 1: 


||o; a ’ A - aG - XaMG\\ Y < K & \a\X 2 + Jiia 2 A||AA4G|| x . (39) 

In fact, if we proceed as in (J5|, Section 2), this result can be improved as follows: there 
exists K 7 > 0 such that, for |a| < 1, 

||a; a ’ A -aG x \\ Y < K 7 a 2 X , (40) 

where 

f N y/l — A 2 f 1 + A 2 1 — A rD 2 

G\\x) = ———exp(^- —x 1 - —x 2 ), x G R . 

Remark that (.C + XM)G\ = 0 and f R2 G\dx = 1. Since G\ = G + XM.G + 0( A 2 ), we 
see that m is compatible with PH). 
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4) Smoothness in x. Standard elliptic estimates imply that u;“’ A (x) is a smooth function 
of x G R 2 for any (a, A) G R x [0, A 0 ]. Indeed, since cv a,x is a stationary solution of ©, 
we have 

= (£ + AA^)- 1 V- (u“’ A u; a ’ A ) . 

It is not difficult to prove that the linear operator (£ + AA4) _1 V is regularizing, hence a 
bootstrap argument shows that iv a ’ x G H k ( R 2 ) for all k G N. One can also prove that all 
derivatives decay rapidly at infinity, so that cu a,x G <S(R 2 ). 

5) Positivity. It follows from the parabolic Maximum Principle JT^j that cn"’ A (x) > 0 for 
all x G R 2 if a > 0. (Similarly, u a,x (x ) < 0 if a < 0, and we already know from (1^ that 
u a,x = 0 if a — 0.) Indeed, arguing as in the symmetric case A = 0, it is not difficult to 
show that, for any initial data u>o G L 2 (R 2 ) nC°(R 2 ), Eq.© has a unique global solution 
u G C°([0, +oo), L 1 (R 2 )) which satishes f R2 u(x, t ) dx = f R2 coo dx for all t > 0. If u>o > 0 
and uq is not identically zero, the Maximum Principle implies that u(x,t) > 0 for all 
x G R 2 and all t > 0, see e.g. (j7j, Section 2.3). Conversely, if a; 0 has non-constant 
sign, the norm ||u;(-,t)|| L i is strictly decreasing in time (for t > 0 sufficiently small), see 
(0, Section 3.1). Since u;“’ A G L\ R 2 ) n C°(R 2 ) is a stationary solution of © satisfying 
J r2 o;“ ,a dx = a, the properties above imply that u> a,x (x) > 0 for all x G R 2 if a > 0. 

5 Stability of asymmetric vortices 

In this final section, we show that the asymmetric Burgers vortex u; Q,A constructed in 
Section Elis a stable solution of © with respect to perturbations in X , provided A > 0 
is sufficiently small. Fix a G R, A G [0, Ao], and consider solutions of © of the form 
u = cd a ’ x + d), u = u“ ,A + u, where u is the velocity held obtained from u by the Biot- 
Savart law ©. Then u satishes the equation 

d t u + u q ’ a • VcD + fi • Vw°’ A + fi • Vch = (C + \M)u . (41) 

If we further decompose a; a,A = aG + w a,x , u“ ,A = av G + v a,A , this equation becomes 

d t Cj + u • Vu> = + XM — aA — A/’ Q,A jd) , (42) 

where A is defined in m and Af a,x is the integro-differential operator dehned by 

A r a ’ x £u = v“’ A • VA + a • Vw a ’ x . 

It is easy to show that the Cauchy problem for ®) or (gH) is locally well-posed in the 
space X. In the symmetric case A = 0, this is proved in py using a larger function space 
(with polynomial instead of Gaussian weight), and the same arguments apply here with 
straightforward modifications. Our goal is to control the behavior of the solutions of (E21) 
in a neighborhood of the origin. An energy estimate yields the following result: 

Proposition 5.1 There exist positive constants K 8 ,K 9 (independent of a and X) such 
that, for any S G (0,1), any solution solution of (El in X satisfies 

4||i.(i)E < -(l-y|l<i(f)|| 2 Y 

whenever Ii 8 X + AVg||a>(t)||yy < S- 
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Proof: From (TI21) we obtain 


— ||w||.Y = (&, (£ + AM - aA - Af a,x )ivJx ~ (^ u ' • (43) 

To simplify the subsequent expressions, we define / = G _1 / 2 d), so that ||d) \\x = \\f\\», 
and we introduce the quadratic form 

£(/) = ./J v/|2+ '€ f ) dx - 

Obviously ||/||| 2 < CE(f), and we have ||/|| 2 4 < Cl|V/||^||a;/||^ 2 < CE(f ), see 0. 
Proceeding as in m we obtain for any 5 G (0,1): 

(a,Cu) x = -U,Lf )y < -4=^11/111= - |b(/) . 

(The bound (fT71) was the particular case 5 = 1/2.) Next, since A4 = \{x\d\ — ^ 2 ^ 2 ) and 
G -1//2 V(D = V/ - f /, we find 

(w,^ < |MU||Vo>|U = ||x/|Ua ||V/ - f/|| L 2 < CTd(/) . 

Moreover, as was observed in the proof of Lemma 12.31 the operator A is skew-symmetric 
in X, hence (a), A Cj)x = 0. 

We now bound the nonlinear term in (TEA . Integrating by parts and using the relation 
VtW 1 = § G -1 , we obtain 

(d), u • VCo)x = / G _1 d)(u ■ Vd)) dx = — - f G^(x ■ u)d) 2 dx . 

J R2 4 J R 2 

Applying Holder’s inequality, we find 

|(d),u-Vu>) x | < 7 / |u||x/|I/I dx < i||u|| L 4||/|| L 4||x /|| L2 < C'||£|| x £(/) , 

4 7r2 4 

where we have used the bound ||u||x,4 < C||d)|| L 4/3 < C||u)||x, see (OHl) and (J5TJ). 

It remains to bound (d),A A u j) x = (d), v q,a ■ Vd;)* + (a;, u ■ Vw a ’ A ) Y . The first term 
in this sum can be estimated in the same way as the nonlinear term above, namely 
|(d), V q,A ■ Vd))* | < C\\w a ' x \\xE(f). For the second term, we argue differently: 

|(d>,u- Vw°’ A )*| < / G _1 |d)||u||dx 

Jr . 2 

< ||V W q ’ a |U||u|| L 4||/|| L 4 < C\\w a,x \\yE(f) , 

where we have used ||u|| L 4 < C||d)|| L 4/3 < C '||/|| i 2 < CE(f) 1 A 2 . Thus, using Proposi¬ 
tion 021 to bound w a ’ x : we conclude that |(d>,A/” Q! ’ A d>)jY| < C\\w a,x \\YE(f) < CAE(f). 
Summarizing, we have shown that 

—-II^H^ + 2 + -^IMU — S^E(G 1/2 d)) , 
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for some constants K$, Ii 9 independent of a, A and of the solution to of (TH1) . This concludes 
the proof. □ 


Theorem 11.21 is now a direct consequence of Proposition 15.11 Given /i G (0,1/2), it 
suffices to choose Ai G (0, A 0 ] such that AT 8 Ai < 8/2 and £ > 0 such that K 9 e < 8/2, 
where 8 — 1 « 2/x. If Co 0 G A" satishes ||u)o||x: < £ and if Cj G G°([0, T*), X) denotes the 
(maximal) solution of (H21) with initial data Co 0 , we dehne 


T = 


supjf G (0,T*) A' 9 ||o;(s)||x < 8/2 for 0 < s < f j G (0,T*] 


If T < Too then Proposition 15.11 implies that ||u)(f)||x < e _/ ^||u)o||.Y for t G [0, T\, hence 
AGg||d>(t)|| jy' < K 9 £ < 8/2 for t G [0, T\, which contradicts the definition of T. Thus, we 
must have T = T* = Too. This means that the solution to is globally defined, and by 
Proposition 15.II ||d;(t)||x < e _/it ||u)o||x: for all t > 0. □ 

We conclude with a few remarks on the basin of attraction of the asymmetric vortex 
and the decay rate in time of the perturbations: 

a) Proposition 15.11 shows that the decay rate in time of pertubations of the asymmetric 
vortex to a,x is bounded from below by [i = |(1 — K 8 A) — \ — 0( A), uniformly in a G R. 
This is consistent with the information we have on the spectrum of the linearized operator 

£°’ a = £ T \M -aA- Af a ’ x , 


acting on the space A". Indeed, differentiating the identity (£ T AA1 )o;“’ a = u“ ,A ■ Vu; a,A 
with respect to x\ and X 2 we obtain 

C a ' x (dno a ' x ) = -^(<9i lu q ’ x ) , jC a ’ x (d 2 io a ’ x ) = -^(d 2 uo a ’ x ) . 

In particular, since d 2 u a,x G X , we see that —i s always an eigenvalue of £ Q,A , hence 
/r < -h-A^ Numerical calculations by Prochazka and Pullin [2] seem to indicate that — 
is always the largest eigenvalue of £“ ,A in A" (for any A). If this was true our arguments 
could be extended to prove existence and stability of asymmetric Burgers vortices for 
all A G [0,1). Put another way, our current limitation on the range of the asymmetry 
parameter is only due to the fact that we do not know how to control the eigenvalues of 
C°' x (except of course for small A). 

b) A remarkable feature of our stability result (Theorem 11 . 21 ) is that it holds uniformly 
for all a G R. In particular, this implies a uniform upper bound on the eigenvalues of the 
linearized operator C a,x . This is definitely compatible with the numerical observations of 
Prochazka and Pullin CHI, but these calculations suggest that our result is perhaps not 
optimal for large Reynolds numbers. According to PI we expect that the eigenvalues that 
are not frozen by symmetries have a real part that converges to — oo as |a| —» oo, which 
could imply a faster decay rate (i and a larger basin of attraction £ for large Reynolds 
numbers. A mathematical understanding of these numerical observations is still lacking. 

c) We chose to consider perturbations Co in the weighted space X because of the “miracu¬ 
lous” fact that the operator A is skew-symmetric in that space. This is why Proposition l5.ll 
holds uniformly for all a G R. However the space X is relatively small since its elements 
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are forced to decay rapidly at infinity in space. Extending the methods developed in [7j 
for the symmetric case A = 0, it is not difficult to show that the asymmetric Burgers vor¬ 
tices are also stable with respect to perturbations in weighted L 2 spaces with polynomial 
(instead of Gaussian) weight. The decay rate in time of the perturbations is still uniform 
in a, but the size of the basin of attraction is a priori not. 

d) In the symmetric case A = 0 it is shown in jj that the Burgers vortex aG is the unique 
stationary solution of © such that u € L^R 2 ) and f R2 u dx = a. Moreover, any solution 
u) G C°([0, +oo), L^R 2 )) of © such that f R2 cv(x,t) dx = a converges to aG in L^R 2 ) 
as t —> +oo. We do not know if such global results can be extended to the nonsymmetric 
case A > 0. 

Acknowledgements. A part of this work was completed when CEW was a visitor at 
Institut Fourier, Universite de Grenoble I, whose hospitality is gratefully acknowledged. 
The research of CEW is supported in part by the NSF through grant DMS-0405724, and 
the work of ThG is supported by the ACI “Structure and dynamics of nonlinear waves” 
of the French Ministry of Research. 


References 

[1] J.M. Burgers. A mathematical model illustrating the theory of turbulence. Adv. 
Appl. Mech., 1:171-199, 1948. 

[2] L. Caffarelli, R. Kohn, and L. Nirenberg. First order interpolation inequalities with 
weights. Compositio Math., 53(3):259 275, 1984. 

[3] E. A. Coddington and N. Levinson. Theory of ordinary differential equations. 
McGraw-Hill Book Company, Inc., New York-Toronto-London, 1955. 

[41 D. G. Crowdy. A note on the linear stability of Burgers vortex. Stud. Appl. Math., 
100(2): 107—126, 1998. 

[5] C. Eloy and S. Le Dizes. Three-dimensional instability of Burgers and Lamb-Oseen 
vortices in a strain held. J. Fluid Mech., 378:145-166, 1999. 

[6] Th. Gallay and C. E. Wayne. Invariant manifolds and the long-time asymptotics 
of the Navier-Stokes and vorticity equations on M 2 . Arch. Ration. Mech. Anal., 
163(3):209-258, 2002. 

[7] Th. Gallay and C. E. Wayne. Global stability of vortex solutions of the two- 
dimensional Navier-Stokes equation. Comm. Math. Phys., 255(1):97—129, 2005. 

[8] Th. Gallay and C. E. Wayne. Three-dimensional stability of Burgers vortices : the 
case of low Reynolds number. Preprint, 2005. 

[9] J. Jimenez, H. K. Moffatt, and C. Vasco. The structure of the vortices in freely 
decaying two-dimensional turbulence. J. Fluid Mech., 313:209-222, 1996. 


16 


[10] E. H. Lieb and M. Loss. Analysis, volume 14 of Graduate Studies in Mathematics. 
American Mathematical Society, Providence, RI, 1997. 

[11] H. K. Moffatt, S. Kida, and K. Ohkitani. Stretched vortices—the sinews of turbu¬ 
lence; large-Reynolds-number asymptotics. J. Fluid Mech., 259:241-264, 1994. 

[12] L. Nirenberg. On elliptic partial differential equations. Ann. Scuola Norm. Sup. Pisa 
(3'), 13:115-162, 1959. 

[13] A. Prochazka and D. I. Pullin. On the two-dimensional stability of the axisymmetric 
Burgers vortex. Phys. Fluids, 7(7): 1788-1790, 1995. 

[14] A. Prochazka and D. I. Pullin. Structure and stability of non-symmetric Burgers 
vortices. J. Fluid Mech., 363:199-228, 1998. 

[15] M. H. Protter and H. F. Weinberger. Maximum principles in differential equations. 
Prentice-Hall Inc., Englewood Cliffs, N.J., 1967. 

[16] M. Reed and B. Simon. Methods of modern mathematical physics. I. Academic Press 
Inc., New York, second edition, 1980. Functional analysis. 

[17] A. C. Robinson and P. G. Saffman. Stability and structure of stretched vortices. 
Stud. Appl. Math., 70(2):163-181, 1984. 

[18] M. Rossi and S. Le Dizes. Three-dimensional temporal spectrum of stretched vortices. 
Phys. Rev. Lett., 78:2567-2569, 1997. 

[19] P. J. Schmid and M. Rossi. Three-dimensional stability of a Burgers vortex. J. Fluid 
Mech., 500:103-112, 2004. 

[20] E. Stein. Singular integrals and differentiability properties of functions. Princeton 
LIniversity Press, Princeton, N.J., 1970. Princeton Mathematical Series, No. 30. 

[21] G. 1. Taylor. Production and dissipation of vorticity in a turbulent fluid. Proc. R. 
Soc. Lond. A, 164:15-23, 1938. 


17 



